Jodar et al. Nanoscale Research Letters 201 3, 8:258 
http://www.nanoscalereslett.eom/content/8/1/258 



0 Nanoscale Research Letters 

a SpringerOpen Journal 



NANO EXPRESS OpenAccess 



Electronic and transport properties in circular 
graphene structures with a pentagonal 
disclination 

Esther Jodar 1 *, Antonio Perez-Garrido 1 and Fernando Rojas 1 ' 2 



Abstract 

We investigate the electronic and transport properties of circular graphene structures (quantum dots) that include a 
pentagonal defect. In our calculations, we employ a tight-binding model determining total and local density of states, 
transmission function and participation number. For the closed structure, we observe that the effect of the defect is 
concentrated mainly on energies near to zero, which is characteristic of edge states in graphene. The density of states 
and transmission functions for small energies show several peaks associated with the presence of quasi-bound states 
generated by the defect and localized edge states produced by both the circular boundaries of the finite lattice and 
induced by the presence of the pentagonal defect. These results have been checked by calculating the participation 
number, which is obtained from the eigenstates. We observe changes in the available quasi-bound states due to the 
defect and the creation of new peaks in the transmission function. 
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Background 

The theoretical and experimental study of properties of 
graphene has attracted the attention of many authors in 
the last few years since a method to isolate single graphene 
layers was developed (the authors Geim and Novoselov 
were awarded with the Nobel prize). These graphene 
sheets may be stable enough to be freely suspended [1], 
which allows us to use them in solid state experiments. 
Besides, the electronic properties of graphene are sur- 
prising: one finds new quasi-particles described by the 
Dirac equation at low energies that behave like mass- 
less particles. This opens the possibility to study quantum 
electrodynamics properties in solid-state devices and to 
carry out new developments, e. g., biosensors (see other 
studies [2-10]). 

The influence of defects and edges in graphene prop- 
erties has been widely studied [11-13]. Other authors 
made similar studies to ours but considered different 
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geometries: Zhang et al. [14,15] worked on transport with 
narrow ballistic ribbon of graphene with zigzag edges 
including topological defects. Carpio et al. [12] studied 
the electronic properties in a similar geometry but with 
dislocations consisting of heptagon-pentagon pairs in an 
hexagon lattice. The aim of this article is to study the 
combined effects, on the electronic properties, of the geo- 
metric structure and the presence of topological defects. 
We then consider a model with a pentagonal defect (discli- 
nation), henceforth PD, at the centre of a graphene sheet 
with a circular shape (see Figure 1). We characterize the 
electronic and transport properties with the local and 
total density of states, participation number and trans- 
mission function. This work can be useful for the search 
of structures suitable for confinement of Dirac electrons, 
which are the basis for the construction of nanoelectronic 
devices with graphene. 

Methods 

Our geometry consists of a finite circular graphene quan- 
tum dot with 1,011 carbon atoms. For electronic trans- 
port, the quantum dot is connected to two semi-infinite 
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Figure 1 Graphene sheet with the topological defect. Schematic 
geometry of the graphene sheet studied in this work. Note the 
pentagonal defect placed at its centre (in red colour). This structure is 
connected to two semi-infinite graphene leads, which are partially 
shown in the figure (red colour). 



Another useful property is the local density of states: 



leads. In Figure 1, we show the quantum dot and, par- 
tially, the semi-infinite leads. We employ a tight-binding 
model that only takes into account one n -orbital per atom. 
The overlap energy between nearest neighbours is taken 
as ^ = 2.66 eV, where second-neighbour interactions 
are neglected. The advantage of using a single-band tt- 
orbital model resides in its simplicity, being the general 
features of electronic transport in very good agreement 
with those obtained by more sophisticated approaches. 
The hamiltonian can then be written as 



H=-tJ2 (£cj + hc), 



(1) 



<lj>cr 



where c\ /q are the creation/annihilation operators of an 
electron in site i. We expand the wave function in terms of 
the site base. = E/^f \i), where a k is the amplitude 
probability that the electron is to be in site i for the eigen- 
state k. We need to solve H\W k ) = E k \^ k ). Four quantities 
are calculated to characterize the nature of the electronic 
and transport properties on two-circled structures, with 
PD and defect-free (ND) structures: the total density of 
states N(E), the local density of states p (i,E), the par- 
ticipation number P(E) and the transmission function 
T{E). 

Electronic properties for the closed system 

The density of states is determined from the energy spec- 
trum as 



N(E) 



1 D 



E k ). 



(2) 



1 ° 

p(i,E) = -J2\(i\*k)\ 2 8(E-E k ), 



(3) 



which measures how each site i contributes to the com- 
plete spectrum. For a fixed £, it characterizes the spatial 
nature of the state: it is localized when only few sites 
contribute to that energy, or extended when more sites 
participate. Finally, the participation number is defined 
as [16] 



P(E) = 



1 



EtiK^lW 



(4) 



It assesses the wave function spreading so it can 
help to find out the localized or extended nature of 
an electronic state. For a completely localized wave 
function ty k (i) is approximately 8^ — > P & 1 while for 
a typical delocalized wave function on D atoms, 
ty k (i) is approximately and then P ^ D. 

Electronic transport properties: Green's function method 

We calculate transport properties for our graphene 
structures such as the density of states and the trans- 
mission function using Greens function method. In 
order to obtain Greens function, we use the following 
expression [17]: 



G(E) = (EI-H-± L -± R ) ) V 



(5) 



k=l 



where El an( ^ Er are the self-energy terms of left and 
right leads, respectively, and H is the Hamiltonian of the 
conductor, i.e., in our case, the circular graphene sheet 
plus a few unit cells of the leads. In our approach, the con- 
tact leads at opposite sides of the circular graphene sheet 
is the graphene sheet itself extended to make the leads 
semi-infinite. This is equivalent to have reflectionless con- 
tacts in macroscopic conductors. Self-energy terms are 
calculated using the prescription Er/lO* j) = ^ 2 Gr/l(1<, 1), 
where Gr/l (k, 1) is Greens function of the semi-infinite 
lead (right or left) evaluated on sites k and 1, which are in 
contact with sites i and j in the circular graphene sheet. 

We only need to calculate Gr/l in the sites in con- 
tact with the conductor. To do that, we use the formal- 
ism developed by Lopez Sancho et al. [18]. This method 
has the advantage that the number of iterations close 
to singularities is very low compared to other trans- 
fer matrix methods, so it converges very fast and has 
been applied to graphene layers by other authors (see 
e.g. [19]). In this scheme, Greens function is Gr/l = 

(eJ — Hoo — Hoi f\ , where Hqq is the Hamiltonian of 

one isolated graphene cell in the lead, and Hq\ is the 
matrix that takes into account the interaction between 
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two consecutive cells. For the calculation of T, we use the 
iterative method described in [18]. 

From Greens function of the graphene structure, we cal- 
culate the transmission function and the density of states 
as [17] 

r(£) = Tr[f L/R G R f R/L G A ] (6) 

N(£) = -^Tr(lmG). (7) 

In Equation 6, G R / A are the retarded and advanced 
Greens functions, respectively, and i\/r = i[ X!l/r — 

-A 

2Jl/r]. We denote the trace of the matrix considered by 
"Tr", which is extended over the whole matrix. 

Results and discussion 

We have obtained different properties of graphene struc- 
tures with and without pentagonal defects, in order to 
evaluate the influence of the defect and the geometry on 
their electronic properties. For the closed structure, we 
have calculated the total density of states, which is shown 
in Figure 2, for both the defect-free structure (dashed 
line) and with PD (continuous line). We see that the den- 
sity for the structure with PD shows a shoulder near 
E = 0, indicating the existence of additional edge states 
induced by the presence of the PD and the circular shape 
of the structure. The behaviour of the participation num- 
ber confirmes these findings (see Figure 3a for the ND and 
Figure 3b for the PD structures). One can observe that 
Ppd < Pnd> which suggests that the PD induces local- 
ized states beyond the known edge states in the structure 
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Figure 2 Density of states for the closed structure. Total density of 
states for the structure with no defects (dashed line) and with a 
pentagonal defect (continuous line). Note the non-null density near 
zero as a manifestation of the edge defects. 
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Figure 3 Participation number for the closed structure. 

Participation number P(E) of the available energy states for the 
structure with no defects (a) and with the pentagonal defect in the 
centre (b). 



with ND. In addition, near E = 0, the structure of the par- 
ticipation number is slightly different in both structures. 
This shows that the nature and the number of the edge 
states are different, which can be seen in the local density 
function for this energy. 

The edge states are localized so only few states con- 
tribute to a certain site; this is shown in Figure 4 for the 
local density of states at E = 0 and p (i, E = 0) for both 
the ND (Figure 4a) and PD (Figure 4b). Clearly, these are 
edge states, and the PD structure shows contribution from 
two zones, compared to the ND structure with one. The 
effect of the PD on the density of states near E = 0 is of 
geometrical nature; the whole structure is affected by the 
presence of the pentagon since it changes the relative ori- 
entation of the edge sites and induces the creation of edge 
states. This has to do mainly with the atom rearrangement 
in the lower part of the structure, which creates new edge 
states and, clearly, the PD sites do not have an explicit 
contribution to such sites. For larger values of E, in the 
local density of p(i,E = 2.6), more sites contribute to that 
energy (see Figure 5). Specifically, we see the contribution 
of sites around the PD as it can be seen in Figure 5b, where 
a star shape appears. The rest of the sites contribute more 
or less similarly to the structure with ND (Figure 5a). 

As the change in behaviour with the presence of PD is 
near zero energy (around the Fermi energy), we concen- 
trate in the analysis of the transport properties around 
such energy. We have also checked our previous results 
in the open structure calculating the density of states 
(Figure 6) and the transmission function (Figure 7). The 
density of states shows several peaks associated with both 
the presence of quasi-bound states (due to the circular 
confinement and the defect) and localized edge states due 
to circular boundaries of the finite lattice. These results 
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Figure 4 Local density of states for E = 0. Spatial distribution of the local density for p(i,E) for the energy E close to zero in (a) a structure with no 
defect and (b) one with the pentagonal defect in the centre. Due to single-bond atoms (see Figure 1 ), the quantum dot is not fully symmetric 
around a central vertical axis. 




Figure 5 Local density of states for E = 2.6 eV. Same as Figure 6 but for the energy E = 2.6 eV. 
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are clearly observed in the peak structure of the transmis- 
sion function (Figure 7), where we observe changes in the 
quasi-bound states available to transport and the creation 
of new peaks in the transmission function. 

Conclusion 

We have investigated the electronic and transport prop- 
erties of circular graphene layers with a pentagonal discli- 
nation. In particular, using a tight-binding model, we have 
calculated the density of states, transmission function, 



participation number and local density of states of the 
structure with and without defects. The density of states 
for the structure with the PD shows several peaks that 
are associated with new localized states, which have been 
checked by calculating the local density of states and the 
participation number. We observe changes in the avail- 
able quasi-bound states due to the defect and new peaks 
of the transmission function. Comparing these results, 
we conclude that there are more quasi-bound states in 
the structure with the defect, states associated with both 




Figure 7 Transmission function for the open structure. Transmission function with the pentagon at its centre (red line) and without it (black 
line). Note the additional peaks for the case with the defect. 
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the presence of quasi-bound states related to the atoms 
belonging to the defect and others due to the circular con- 
finement and edge states due to circular boundaries of the 
finite lattice and the defect. 
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